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ABSTRACT 

A hadronic chiral SU(3) model is applied to neutron and proto-neutron stars, 
• taking into account trapped neutrinos, finite temperature and entropy. The tran- 

I ■ sition to the chirally restored phase is studied and global properties of the stars 

;-H ■ like minimum and maximum masses and radii are calculated for different cases. 

^ ■ In addition, the effects of rotation on neutron star masses are included and the 

conservation of baryon number and angular momentum determine the maximum 
■ frequencies of rotation during the cooling. 

> 

2^ ■ Subject headings: neutron stars, proto-neutron stars, chiral restoration, star cool- 

ing, star rotation 
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O ■ 1. Introduction 

^ . After a core-collapse supernova explosion of a star with a mass smaller than about 20 

. solar masses, the remaining star, initially called proto-neutron star, is left with a very high 
temperature of up to 50 MeV. Due to their short mean free path especially electron neutrinos 
stay trapped for roughly 10 to 20 seconds and develop a finite chemical potential. In this 
stage the structure of the star can be divided into a core region that will be modelled in this 
article, and an envelope with higher entropy. Maximum and minimum star masses as well 
as rotational constraints on them will be considered. Finally, the corresponding properties 
of the cooled neutron star will be determined as well. 

As a direct solution of the QCD equations for a high-density hadronic environment is 
currently out of reach, the core region of the star will be described using an effective SU(3) 
chiral model including hyperonic degrees of freedom. The effect of strange hadrons on the 
properties of neutron stars has been discussed in the past by a number of authors including 
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the investigation of kaon condensation (IKaplan et al.l (19861): 



of th e effect of hyperons in neutron s tars (Glendenning et al 



^ons et al. ( 


200G)) and studies 


( 


1991 




2001 


) 


Jha et al. ( 


2006. 



20071 ) ) and references in iPage et al.l (120061 ) ) . Here in a general approach the baryons and 
mesonic fields are introduced as fiavor-SU(3) multiplets. The baryons obtain their masses 
through their coupling to the scalar meson fields via spontaneous symmetry breaking. The 
parameters of the theory are fixed to reproduce the hadronic vacuum masses as well as the 
properties of ground state nuclear matter. The model is extended to include effects from 
the lowest spin-3/2 baryonic multiplet. The r nodel has been successfu l ly app lied to the de- 
scription of nuclear matter and finite nuclei (iPapazoglou et al. Jl998l . Il999h as well as to 
the high-temperature an d high-density environment of an ultrarelativistic he avy-ion collision 
(IZschiesche et al.l (|2007 ) ). A s tudy of neutron stars had also been performed (IHanauske et al. 



( 119991 ): ISchramm et al.l (l2003l )). In the following these calculations are extended to proto- 
neutron star environments with finite temperature and neutrino chemical potential, addi- 
tionally considering the population of higher baryonic resonances. In a further refitting of the 
basic model parameters, a study of the maximally achievable star masses in this approach 
is presented, while still maintaining a phenomenologically reasonable description of ground 
state nuclear matter. 

The article is organized as follows. In section 2 the hadronic degrees of freedom and 
the model Lagrangian will be briefly presented. In section 3 the model will be applied 
to zero-temperature neutron stars with special emphasis on the effect of non-linear meson 
interactions. The properties of proto-neutron stars at finite temperature, entropy and lep- 
ton number are shown and discussed in section 4. Finally, the effects of rotation and its 
consequences on the cooling of the star are analyzed and conclusions are presented. 



The Hadronic Model 



Since in a high-density environment like the one present in the center of neutron stars 
the baryonic chemical potential is high enough to create particles beyond the lowest SU(2) 
multiplet of the nucleons, which comprise hyperons (A, S, S) and possibly other resonance 
states (A, S*, H*, Q), the hadronic model is based on a fiavor-SU(3) description including 
hyperons and strange mesons as part of the basic hadronic SU('3) r nultip le ts, which were al so 
partially taken into account in previous works (ISchaffner et al.l (119961 ): |Jha et al.l (120071 )). 
The baryons interact via exchange of scalar {a, S, (, x) cind vector mesons (p, uj, 0). The 
isovector vector (p) and scalar (5) mesons are included in a natural way, as an essential 
ingredient for reproducing the phenomenological value for the asymmetry energy of nuclear 
matter. A scalar isoscalar (dilaton) field (x) acts as an effective gluon condensate. On the 
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leptonic side, electrons and muons assure charge neutrality and neutrinos are considered in 
the case of proto-neutron stars. 

The Lagrangian density of the chiral model in mean field approximation used in our 
calculations reads: 

LmFT = Lxin + Lsscal + Lsvec + Lgcal + L^ec + LgB-, (1) 

where besides the kinetic energy term for baryons and leptons, the terms: 

Lsscal + Lbvcc = -^i^i [giuj^o^^ + 5'i07o0 + giploT3p + m*]%l)i, (2) 



(3) 



(4) 



'SB 



(5) 



represent the interactions between baryons and scalar mesons and between baryons and 
vector mesons, the self interactions of scalar and vector mesons and an explicitly chiral 
symmetry breaking term, responsible for producing the masses of the pseudo-scalar mesons. 
In Eqs. (2-5) the mesons are treated as classical fields within the mean-field approximation. 
Only the O*'* component of the vector meson survives, i.e. u = Uq, p = p^, cf) = 0o- and ( are 
the scalar fields corresponding to the non-strange and strang e quark condensate , respectively . 
A detailed discussion of this Langrangian can be found in iPapazoglou et al.l (Il998l . Il999l ). 
The baryon masses are generated by the scalar fields except for a small explicit mass term 
equal to 6m ~ 150 MeV for nucleons. The effective masses decrease at high densities with 
decreasing scalar fields as the chiral symmetry is partially restored. At low densities they 
reproduce the experimentally known baryon masses: 



m 



(6) 



The coupling constants used to calculate proto-neutron star properties {gNuj = 11.9, 
gN<p = 0, gNp = 4.03, gNa = -9.83, gm = -2.34, g^c = 1.22, ko = 2.37, h = 1.40, fca = 
—5.55, fcs = —2.65, A;4 = —0.23, e = 0.06, g^ = 38.9) are chosen to reproduce the vacuum 
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masses of the baryons and mesons, the nuclear saturation properties (density pB = 0.15fm~ , 
binding energy per nucleon B/A = —16.00 MeV, compressibihty K = 297.32 MeV), the 
asymmetry energy {Egym = 32.50 MeV), and reasonable values for the hyperon potentials 
([/a = -29.41 MeV, f/s = 20.39 MeV, = -10.09 MeV). The vacuum expectation values 
of the scalar mesons are constrained by reproducing the pion and kaon decay constants. 



3. Neutron Stars 

In the calculation of zero-temperature neutron stars the influence of the structure of the 
self interaction term of the vector mesons is investigated , as in dense systems baryonic vector 
densities, and therewith the mean fields of the vector mesons, become especially important 
for the equation of state of hadronic matter. The fourth-order self-interaction term Lt,ec,4 of 
the vector mesons can be written in different forms in a SU(3)-invariant way. To study the 
difference in the result three separate coupling schemes are considered, 

UecA = -9^[Tr{V)]y4 (a), - g,[Tr{V')]' (6), - g,2Tr{V') (c) , (7) 

where V stands for the 3x3 matrix of the vector meson multiplet, which reduces to a diagonal 
form in the mean field limit, i. e. V" = {{u + p)/V2, {cu - p)/V2,(f)). The mass-radius 
relation of the respective neutron star for the cases (a) to (c) was calculate d by solv i ng th e 



Tolman-Oppenheina e r-Vol kov (TOV) equations for static spherical stars (ITolmanl (119391 ): 



Oppenheimer et al.l (119391 )). The results are shown in Fig. [T]in the SU(2) limit of proton, 
neutron and electron matter. A strong coupling of the u and p meson present in the cases 
(b) and (c) leads to smaller star masses. In the following calculations the non-linear coupling 
(a), which does not generate a. p — u coupling, is used. This allows for more massive neutron 
stars and is also in general agreement with the observed small mixing of the two mesons. 



The modification of the original parameter set CI of the model, used in lPapazoglou et al. 



(119991 ) was done in order to investigate the maximum neutron star masses th at can be 



achieved in the model (similar studies in a different approach have been done in iJha et al. 



(120071 )). while still reproducing the hadronic masses in the vacuum as well as reproducing 
the phenomenological values of the basic nuclear matter ground state properties as listed 
in the previous section. Within the model different situations can be analyzed by including 
the whole baryon octet or, in addition, the baryon decuplet. In practice, the only baryons 
present in the star besides the nucleons are the A and S~ in the first case (Fig. [2]) and A, 
^-,0,+,++^ in the second case (Fig. [3]). In the presence of resonances the A~ particle replaces 
the E~, as its effective mass drops faster with density compared to the S~. As can be seen 
from Fig. Hlthe inclusion of new particles, i.e. new degrees of freedom, softens the equation 
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of state (EOS). The same effect would be observed in the symmetric case, when there is 
no net isospin instead of no net electric charge, in the presence of more massive degrees of 
freedom. Although the symmetric case has no relevance for neutron stars due to their high 
intrinsic asymmetry this case is important for example in heavy ion collisions. 

The softening of the equation of state causes a decrease of the respective neutron star 
maximum mass (Fig. [5]). Besides using the high density EOS for the interior of the star, 
in order to calculate the mass-radius diagram, an i nner crust, an outer crust and an atmo- 



sphere have been considered following iBaym et al.l (jl97ll ). The studies of the proto- neutron 



star properties are restricted to the model including the lowest multiplets, i.e. including 
the baryon octet, thus avoiding the uncertainties related to the largely unknown coupling 
strengths of the baryonic decuplet. In this case, it is possible to describe stars with masses 
higher than M = 2Mq and still take into account heavy baryonic degrees of freedom. 

The transition to the chirally restored phase for any of the considered self-interactions 
and sets of baryonic degrees of freedom turns out to be a cross over. This happens due to 
the requirements of beta equilibrium and charge neutrality that make the different isospin 
states of baryons with the same vacuum mass appear at different densities of the star, 
thus smoothing out their effect. This transition can be seen in the behavior of the scalar 
condensate a used here as the order parameter for the transition. In Fig. [6] the condensate 
was plotted against the star radius showing that in this model the chiral symmetry is partially 
restored in neutron stars. 



4. Proto-Neutron Stars 

Right after the supernova explosion, due to the neutronization of the matter, the star 
contains an abundant number of neutrinos that are trapped in the system. The temperature 
of the star can reach values of up to 50 MeV. In order to determine the matter properties 
under these conditions the thermodynamical potential of the grand canonical ensemble is 
solved. It is defined as: 

£2 T T T T 

y J-'scal J-'vec J-'SB ^vac 

TTE.. ln(l ± e-^(^rW-Mr)), (8) 

where i denotes the fermion type (including leptons), 7^ the fermionic degeneracy, E*{k) = 
a/aJM-^P the energy, and fi* = fJ'i — gu^^o — gpPosTsf^ — g(i)4>o the effective chemical potential 
(a vanishing chemical potential for muon neutrinos has been assumed here). In the case of 
leptons m* = m and /i* = /x. The single particle energy is given by Ei{k) = E*{k) + 
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Qi^uj + gppr^/2 + g^(j) and the entropy per volume per baryon is defined as s = S/V/ ps = 
-dn/dT\v,^,/pB- 

To study this complex system two different features are taken into account separately: 
finite temperature and high lepton number. In a static approximation of the evolution of 
the star three different approaches are considered as they have appeared in the literature: 



constant temperature: in this case the whole star is considered at the same temper- 
ature, which is unrealistic, and the maximum mass of the star is higher for higher 
temperatures because of the thermal effects on the binding part of the mass. In this 
case the entropy per baryon remains constant with the increase of density except for 
small densities as can be seen in Fig. [3 This effect comes from the fact that even when 
the baryon density tends to zero, the electron-positron pairs present at finite temper- 
ature still contribute for the entropy. The transition to the chirally restored phase 
becomes smoother with the increase of temperature since the Fermi surface becomes 
less important (Fig. [8]). 



metric dependent temperature: in this discussed in iGondek et al.l (119971 ) , the 



temperature is defined at an infinite distance from the star and the thus defined tem- 
perature increases as the gravitational field created by the presence of the star mass 
becomes higher: T = Tqo/a/^, where goo is the 00 component of the metric tensor. As 
the gravitational field increases toward the center of the star the density also increases 
as can be seen in Fig. [9] for Too = 15 MeV, but the increase of the temperature from 
the center to the edge of the star is not too pronounced. In this case, the maximum 
mass of the star is higher for higher temperatures Too. 

constant entropy: in this case the star is considered to have a constant entropy through- 
out, which t o _some degree agrees with dynamical simulations of the stellar evolution 



(jStein et al.l (120061 )). The temperature is higher in the center of the star and colder at 
the edge as can be seen in Fig. [10] and the maximum mass of the star is higher for 
higher entropies also because of the thermal effects on the binding part of the mass. 



The trapped neutrinos with a chemical potential id^, are included by fixing the lepton 
number defined as Yi = (pe + Pue)/PB- In consequence there will be a large number of 
neutrinos in the star but also an increased electron density. Therefore, demanding charge 
neutrality, the proton density increases, and with higher proton density, the star becomes 
more isospin symmetric and the neutron Fermi energy decreases. Thus the increase of 
lepton number softens the EOS and consequently the maximum mass of the neutron star 
gets smaller. It is important to keep in mind that this result is dependent on the chosen 
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parameter set and consequently on the particles present in the star. For parameter sets 
that allow the hyperons to appear at lower densities, the high proton density delays their 
population allowing the neutron star to be more massive. 

These two features can be put together to describe the evolution of the star for the 
extreme cases: constant temperature and fast increasing temperature throughout the star. 
After the supernova explosion, the star is still warm so the temperature is fixed to T = 30 
MeV or the entropy per baryon is fixed to s = 2 (the temperature increases from at the 
edge up to 50 MeV in the center as in Fig JTU]) . The star still contains a high abundance of 
neutrinos that were trapped during the explosion so the lepton number is fixed to a typical 
value of Yi = 0.4. After several seconds (10-20) the neutrinos escape and /3-equilibrium is 
established. After about a minute, the temperature of the star has dropped below 1 MeV. 

The balance between the effects of temperature/entropy and lepton number is very 
delicate and depends on the parameters of the model. The first line in table [1] shows the 
maximum star masses for different entropies and different temperatures. As can also been 
seen in Figs. [TT] and [T2] the intermediate step of the evolution with s = 2 and /ij, = is more 
massive than the first one. But this calculation does not take into account that the baryon 
num ber does not incr e ase a s the star gets colder, otherwise it would collapse into a black 



hole (iTakatsuka et al.l (Il995l )). If the baryon number is fixed starting with the values for the 
star with s = 2 or T = 30 MeV, resp., and Yi = 0.4, the stable solutions of the colder cases 
have a smaller mass than the warmer case as can be seen in the second line of Tab. [H 

In stars with finite temperature/entropy a crust of high entropy has been used since 
the shock wave created during the supernova explosion leaves the outer region with a much 
higher entropy than the rest of the star and it remains warmer for longer time serving as 
an insulating blanket, which d elays the star from co ming to a complete thermal equilibrium 



with the interstellar medium (ILattimer et al.l (jl99ll )). In this case the crust is stiff enough 



to generate massive stars for small central energies resulting in big radii. Because of the 
use of a warm crust it is also possible to stipulate a minimum mass for each case (fourth 
line of Tab. [1]). For the assumption of constant temperature the crust used has an entropy 
s = 5 so that the inner and outer EOS could be continuously c onnected but as a resu lt the 



minimum mass found is far from a realistic value (M ~ 1.18M„. iFaulkner et al.l (120041 )). For 

the constant entropy case the crust used has a value of s = 4. 
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Rotation 



As the proto-neutron star cools down, its possible radii and masses change (Tab. [T]) 
which means that its moment of inertia changes and for determining the ma ximum frequency 



of the star one has to take into account angular momentum conservation (ITakatsuka et al. 



( Il995l )). The maximum frequency (Kepler frequency) with which a star can rotate without 
starting to expel matter on the equator, has been determined by including monopole and 
quadrupole corrections to the metri c due to the rotati o n and solving the self-consistency 
equation for Qk as it was derived in iGlendenning et al.l (119941 ) . In the current calculations 
rotational instabilities, that would potentially reduce the maximum rotational frequency, 
are not taken into account. The rotation of the star generates a modification of the metric 
and the higher the rot ational frequency, the higher the mass and radius of the star become 
( jSchramm et al.l (120031 ) ) as can be seen in Fig. [131 The maximum masses calculated in this 



approximation differ by less than 1% from full calculation results ( IStergioulas et al. I (119951 )) 
for frequencies around 716 Hz. For this set of parameters the increase in the maximum 
masses of the stars from u = to u = u^, fixing the baryon number to the value for zero 
frequency, is smaller than 5%, different from the 15% mass increase in the case that the 
baryon number is not fixed. But this situation can be identified as the spin down of a cold 
star with a certain baryon number (in this case 0.28 x 10^^ N) that continues until it emits 
all its energy and stops rotating. In this point the star is considered dead. 

To consider the rotational effects in the cooling from the very first stage of the neutron 
star life, besides the baryon number, the angular momentum has to be taken into account 
since it remains about constant in the first seconds of the cooling. This new constraint lowers 
even more the possible maximum mass and Kepler frequency of the star as can be seen in 
Tab. [H As a consequence of that, some stars are stable during a stage of cooling but not in 
the subsequent one, causing matter to be lost during the process. 



Conclusion 



For the first time an effective model including hyperon degrees of freedom could generate 
a neutron star with a mass M ~ 2.1 in agreement with the most massive one observed 



but still under discussion (INice et al.l (120051 )). The inclusion the baryonic spin-3/2 multiplet 
as additional degrees of freedom still leads to a quite heavy star of about 1.9 solar masses. 
This mass might be lowered when the cooling process of the star is examined and the baryon 
number is fixed in the early stage of the proto-neutron star life. The maximum mass of the 
star increases again when rotation is included to values beyond M = 2.0Mo for a freq uency 
of u = 1122Hz, the highest frequency measured for pulsars so far (iKaaret et al.l (120071 )) but 
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because of the baryon number and angular momentum constraints, a range of masses in fast 
rotating stars leads to unstable states and shedding of matter. 

The transition to a largely chirally restored phase takes place in the core of the neu- 
tron star within the chiral model, although for all scenarios that were considered here, the 
restoration occurs via a smooth cross over due to the requirements of beta equilibrium and 
charge neutrality. 
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Fig. 1. — Star mass versus radius for different vector meson self interactions. 




Fig. 2. — The composition of neutron star matter with hyperons. 
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Fig. 3. — The composition of neutron star matter with hyperons and resonances. 
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Fig. 4. — Binding energy per nucleon versus baryon density. 
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Fig. 5. — Star mass versus radius for different compositions. 




Fig. 6. — Scalar condensate versus star radius. 
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Fig. 7. — Entropy per volume per baryon versus baryon density for different temperatures. 
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Fig. 8. — Scalar condensate versus chemical potential for different temperatures. 
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Fig. 11. — Star mass versus radius for different moments of tfie star evolution defined by 
temperature and lepton number considering the baryon octet, e, /i. 
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Fig. 12. — Star mass versus radius for different moments of the star evolution defined by 
entropy and lepton number considering the baryon octet, e, /i. 
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Fig. 13. — Star mass versus rotational frequency of the star for different compositions. 
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Table 1: Maximum masses with and without fixing the baryon number, radii for the maximum 
masses with fixed baryon number and minimum masses with and without fixing the baryon 
number for different moments of the star evolution. 



stages 


s = 2 


s = 2 


s = 


T = 30 MeV 


T = 30 MeV 


T = MeV 




Yi = 0.4 


P equiL 


P equil. 


Yi = 0.4 


/? equiL 


/? equil. 




2.05 


2.07 


2.05 


2.08 


2.08 


2.05 


fixed Mb 


2.05 


2.01 


1.96 


2.08 


2.03 


1.99 


R{km) 


14.00 


15.87 


12.44 


16.09 


17.23 


12.33 


MminiMo) 


1.07 


1.07 


0.02 


2.02 


2.02 


0.02 


fixed Mb 


1.07 


1.07 


0.94 


2.02 


2.02 


1.77 



Table 2: Maximum frequencies and masses with and without fixing Mb and L for different 
stages of the cooling. 



stages 


s = 2,Yi = 0.4 


s = 2, /3 equil. 


s = 0, j3equil. 


s = 0, (5 equil. 


vk{Hz) 


1452.56 


1424.48 


1515.86 





Mmax{Mo) 


2.31 


2.33 


2.34 


2.06 


fixed Mb 


1228.23 


1110.28 


1162.03 







2.25 


2.19 


2.14 


2.06 


fixed L 


1195.30 


1095.90 


1162.03 






2.24 


2.19 


2.14 





